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We study bosonization of the sine-Gordon theory in the presence of boundary interac- 
tions at the free fermion point. In this way we obtain the boundary S-matrix as a function 
of physical parameters in the boundary sine-Gordon Lagrangian. The boundary S-matrix 
can be matched onto the solution of Ghoshal and Zamolodchikov, thereby relating the 
formal parameters in the latter solution to the physical parameters in the lagrangian. 
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1. Introduction 

It is now understood that for the integrable quantum field theories in two dimen- 
sions, certain boundary interactions which preserve the integrability are possible 0. Apart 
from the theoretical interest, these models have interesting applications for example to ID 
impurity problems and edge excitations in fractional quantum Hall states ||] . 

In pL boundary scattering matrices were computed for the sine-Gordon model. There, 
a general solution to the boundary Yang-Baxter equation and the crossing-unitarity condi- 
tions was obtained which depends on two formal parameters. An important open problem 
in this work was the relation of these formal parameters to the physical parameters in 
the lagrangian. In this paper we obtain this relation at the free fermion point of the 
sine- Gordon theory, and see that even here this relation is somewhat non-trivial0. (In 
applications to the fractional quantum Hall effect, the free fermion point corresponds to 
filling fraction v = 1/2.) Our computation involves understanding bosonization in the 
presence of boundary interactions, which to our knowledge has not been studied before. 



2. Bosonization with Boundary Interactions 

The sine-Gordon theory with an integrable boundary interaction is defined by the 
Euclidean action 

S = -L J dxdt^ (d x $ds$ - 4Acos/3$) - J- J dt cos - </> )/ 2 ) , (2- 1 ) 

where z = (t + ix)/2, ~z = (t — ix)/2. Here, A and g are dimensionful parameters, (3 is 
the coupling and <po is a constant parameteri. The boundary is taken to be the time axis 
— oo < t < oo at x = 0. 

One may view the above action as a boundary perturbation of a conformal field theory. 
Conformal field theories with boundary conditions that preserve conformal symmetry were 
studied by Cardy||. Setting g = A = 0, and requiring the boundary terms to vanish in the 
variation of the resulting conformal action yields the 'free' boundary condition: = 
at x = 0. In the conformal limit, $ = (p(z) + 4>(z), and this implies 

free : <f>(z) = $(z) -a (x = 0), (2.2) 

1 After this work was completed we learned that similar results were obtain in [jjj). 

2 (5 is related to the conventional coupling /3 by (3 = (3/ \/4tt, and (f>o in this paper differs from 
the conventions in by <f>o — » 
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where a is a constant which represents data not specified by the action (|2.1| ). 

Consider now the limit g — > oo, where the bulk is still conformal with A = 0. In 
this limit, vanishing of the boundary terms which arise upon varying the action requires 
sin(/3($ — 0o)/2) = 0. This preserves conformal invariance, and in terms of the chiral 
components of the scalar field implies the 'fixed' boundary condition: 

fixed : <j>{z) = -(p(z) + (f> + (x = 0), (2.3) 

P 

where n is an integer. 

For arbitrary g, the theory interpolates between the above free and fixed boundary 
conditions. This model was further studied in |7j , and when A = studied in || . Interest- 
ingly, when (3 = the model is conformally invariant for all g (when A = 0) . 

It is well-known that at (3 = 1, the bulk theory is equivalent to free massive charged 



fermions[10| [11]. We now describe how to extend this bosonization to include the boundary 
interaction. Consider first the situation where the bulk is massless. Let if)±(z), if}±{z) be 
the left and right chiral components of the Dirac fermion with U(l) charge ±1. The 
bosonisation relations read tp±{z) = exp(±z</>(z)), tp±(z) = exp(=|=i</>(z)).i With these 
bosonisation relations, A is identified with — m, the soliton mass. The relative minus sign 
in the exponent for left versus right comes from the fact that the U(l) charge in the 
massless limit is —i{§ ^j9 z (p — § -^d^ifi). Thus, in terms of the fermions, the fixed and 
free boundary conditions read at x = 0: 

free : ih± = e Tl<J tb^ 

- (2-4) 
fixed : V± = e ± ^° ip± 

One sees that free boundary conditions break U(l) symmetry whereas fixed boundary 
conditions preserve it. 

We formulate the fermionic action as a perturbation of the free boundary condition. 
The action which enforces the free boundary condition on the fermions is the following: 



1 f 1 — — — — 

S'free = I dxdt- (ifj+dfl/j- + i()_d&l>+ + t(f + dgip_ + ^_d z ^ + ) 

- ^ y eft (e iCT V+^+ + e"^-^_) . 



(2.5) 



3 These bosonisation relations depend explicitly on the exact form of the Dirac action. However 
with these conventions, charge conjugation is implemented simply by taking $ — > — <E>. 
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Namely, the boundary terms in the variation of the above action yield ( |2.4| ). 

One can now express the boundary perturbation in terms of the fermion fields. The 
chiral components of the free boson have the following standard mode expansions in the 
conformal limit :§ 

—i(j)(z) = — z~ n — cio log z — 5o 
^ n 

-_ V-«n_ _ _ ^ 

— iMz) — y — z n — «o log z — «o 
n 

with [a n , a m ] = \a n , a m ] = n5 nj _ m , [oto, 5?o] = [«o,«o] = 1- Since we view the in- 
terpolating theory as a perturbation of the free boundary condition, </>, </> must satisfy 
d z (j)(z) = dj(f)(z) at x = 0, which implies a n = a n , leaving no constraints on 5 ,ao- On 
the boundary one has: 



$ = <j)(z) + (j){z) = 20 - i(a - a ) 
= 2(f) + i(a - a ). 



(2.7) 



Therefore at x = 0, 



cos 



(^-^j = (V> +G _ + ^_a+) e~^/ 2 + (a + </>_ + a_^ + ) e^ 2 {2.1 



where a± are zero mode operators: 

o± = - exp(±z(7r/2 + <f>(z) - 0(Z))/2) = - exp (±{m /2 + a Q - & )/2^ . (2.9) 

The additional factors of 7r/2 arise in combining the ^'s and the a's via exp(A) exp(S) = 
exp(l/2[A, B}) exp(A + B). This leads us to consider the complete action 

S = Sfree "T" 'S'mass + 'S'int. (2-10) 

where the interpolating term is 



S'int. = j dt 

1 

~ 2^ 



(^+a_ + e"^ o/2 + (a+^_ + a_V>+) e^ o/2 



dt a+dtd- 



[2.11) 



See for example [12] 
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and S mass is the bulk mass term — J dxdt(i/j-i/j + — vfj_vfj + ). The action ( |2.10| ) is com- 
pletely analagous to the Ising case considered in [|l]. From the action and the dimension 
of the fermion fields one sees that g has dimensions of y / mais. 

Varying the action with respect to the fermion fields and the zero modes a±, the 
boundary terms, upon elimination of the zero modes, yield the following interpolating 
boundary conditions at x = 0: 

+ e^-^tp- - e~ ia i)_ - e l ^ + = 
id t (ip- - e ia ^ + ) - g 2 (^_e"^° - </>-) = (2.12) 
id t (^+ - e~^_) + g 2 0+ - e^> + ) = 0. 

By taking the g — > (g — > oo) limit, we recover the free (fixed) boundary conditions. 



3. Derivation of the Boundary S-Matrix 

Generally, given a theory with particles A\(6), an integrable boundary perturbation 
produces scattering via the boundary S-matrix, R b a {9)\ 

Al(0)B = R b a (e)Ai(-6)B, (3.1) 

where B is the boundary operator. In our case, the theory has two particles, A\(6) 
(solitons) and A_ (9) (anti-solitons) , and the above equation reads 

aUo)b = p+(e)AU-e)B + q+(9)aI(-9)b-, 

(3 2) 

Al(0)B = Q-{6)A\(-6)B + P-(9)AI(-6)B. 

P + (9) and P~ (9) represent soliton-soliton and anti-soliton-anti-soliton scattering whereas 
Q + (9) and Q~{9) represent soliton-anti-soliton scattering. Because the boundary interac- 
tion, in general, is not U(l) preserving, the appearance of Q + (9) and Q~(9) is expected. 
Only when g —>■ oo should Q + (9) and Q~{9) vanish. 

To derive the pieces of the boundary S-matrix, we interpret the boundary conditions 
in ( p.!2| ) to vanish when acting on B: 

^ + e ,( ^V--e" lff i-e*?JB = etc. (3.3) 



4 



We then substitute the following mode expansions for the fermions: 

^ + = ^rJ°° ^L e 9 ' 2 ^A.(6)e-^ ze0+ ^ - A\{9)e mi - ze6+ze ~^) , 



-oo 



,x w (3-4) 
V+ = -iyfri J ^~ 9 ' 2 (A.(e)e-^ ze +ze ) + A\.(6)e m < ze +Te ^) , 



-oo 

t 



where ip\_ = ?/>_ and ip + = ip_. The A's satisfy the following anti-commutation relations: 

= {A_{e),Al{e')) = A^8{e-e'). (3.5) 

The normalization of the mode expansions is fixed by insisting the two point functions (in 
the bulk) satisfy 

(0\ip_(z,z)'ip + \0) = -2zmK (mr), 

(0\^{z,z)^+\0) = 2my-Ki(mr), 

where r/2 = and the K's are standard modified Bessel functions. Performing the 
substitutions, we find 



p+(0) = (cosh(6>) - I cosh(6> + z0 o )) 
p-(0) = (cosh(6>) - I cosh(6> - z0 o )) /D{9), 
Q+(0) = - % -e~ irJ sinh(29)/D{9), 
Q~(0) = - % -e ta smh(29)/D(9), 



(3.7) 



where D(9) and 7 are given by 

7 = #V TO ; 

DW = ncosh (»+**>+**/*) sinh f*-^-W2\ _ cosh * m , (a8) 



We see that Q + (9) and Q (6*) differ by a phase. To remove the phase, we apply a U(l) 
gauge transformation to A + and A_: 

-> e - <<T/2 ^ , Al^e ia/2 Al. (3.9) 

This leaves P+(0) and P~{9) invariant, but takes Q+(9) -> e i<T Q+(#) and Q - (0) -»• 
e- lfJ Q-(^). Then 

Q+(0) = Q~{9) = - % - sinh(20)/Z>(0). (3.10) 



In what follows, we assume this gauge choice has been made. 

If <po = 0, we expect charge conjugation symmetry to be restored, as indeed it is. 
We also see that the R-matrix is independent of the sign of g. Changing the sign of g is 
equivalent to the shift <p$ — > 0o + 27r. But the subsequent shift $ — > $ + 2n restores the 
boundary term and leaves the bulk terms unchanged, leading to the independence from 
the sign of g. 

Taking the appropriate limits, we can find the boundary S-matrix in the free and fixed 
cases. For the free case we find 

free : P ± (9) = -sech(0), Q±(0) = itanh(0). (3.11) 

In this limit, charge conjugation symmetry, $ — > — <E>, is respected. As a result, P + (0) 
equals P~{9). For the fixed case we obtain 

>±ta\ u f 0±i(/>o-in/2\ , , (0 =F ifo + in/2\ . 



fixed: P ± (9) = -cosh ^ ^ ^j/cosh^ ' ^ J> Q (*) = °- 

(3.12) 

As expected the soliton-antisoliton scattering amplitudes vanish in this limit where the 
boundary interaction preserves the U(l) charge. 



4. Analytic Structure of the Boundary S-Matrix 

Like the bulk S-matrix, the poles of the boundary S-matrix, R, provide information 
on bound states. However in the case of R, the situation is more complicated. Poles 
in R can either be indicative of bulk bound states interacting with the boundary or of 
boundary bound states, which are effectively excitations of the ground boundary state. In 
the former case, poles in R will both be found at = in/2 and at in/2 — 9b, where 9b is 
the location of the pole in S corresponding to the bulk bound state. Such poles necessarily 
imply the addition of zero- momentum states to the boundary state \B). In the latter case, 
the poles may appear anywhere in the region < 9 < in/2. If the pole in this case appears 
at 9 = in/2, a zero-momentum state will, as before, be found in \B). Because the bulk 
S-matrix at the free-fermion point has no pole structure, we are only faced with the latter 
situation. 
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The poles in the range < 6 < in/2 appearing in the free-fermion boundary S-matrix 
depend on the parameters 7 and 0o- There are four areas in this parameter space: 

i) ^ cos0 o > 1, cos0 o < 1 : u = sin -1 (7/4 - (7 2 /16 - ^ cos0 o + l) 1 ^ 2 ) ; 



2 



VI 



cos 0o > 1) cos 0o = 1, 7 < 4 : u = sin 1 (7/2 — 1); 



(4.2) 



7 

Hi) —cos 0o > 1, cos 0o = 1, 7 > 4 : no poles; 

7 

iv ) < — cos 0o < 1 : no poles; 

(4.1) 

where u = —id is the location of the poles. In the free and fixed cases, this pole structure 
reduces to 

fixed: u = 7r/2-|0 o |, < |0 o | < tt/2; 
free : u = n/2. 

In all cases the poles are simple, and we are not burdened with interpreting more compli- 
cated resonance phenomena. 

The energy, E a , of these boundary bound states, B a , is given by 

E a — Eq = m cos(u a ), (4-3) 

where m is the soliton mass, Eq is the energy of the boundary ground state, and u a is the 
location of the pole in the boundary S-matrix corresponding to the state. We expect the 
boundary bound state B a to be stable if E a — Eq < m; that is, if there is no chance of the 
boundary emitting a zero-momentum soliton. We thus see all the boundary states to be 
stable. 

When boundary bound states are present, scattering off the boundary ( j3?ll) needs to 
be described by the more general equation 

Ai(6)B a = R b a i(6)Al(-6)B^ (4.4) 

where {B a } are all the possible boundary states. By energy conservation, we expect 
Rip = if E a 7^ Ep. Thus for all but the free case and the fixed case with 0o = 0, the 
boundary states are not coupled. With such independence, we expect = R h ^a as the 
constraints among the free-fermi fields (|2.12[) must annihilate B a for every a. 

The interpretation of the boundary states in these non-degenerate cases is straight- 
forward. Classically, the asymptotic behaviour of the ground state in the bulk is described 
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by $ — ► (2n + 1)7T as x — > — oo (the bulk classical potential is V(3>) = j^cos($); recall 
A = — m and so A is negative). The difference in energy between these states arises from 
a boundary term. Provided 7 is large enough, the boundary term forces $ ~ cpo at x = 0. 
If 4>o is in the range < <po < n/2, the boundary bound state can be viewed as the state 
where $'s asymptotic value is —n (for the ground state the asymptotic value of $ is n). 
If — n/2 < <fi < the situation is reversed: the ground state corresponds to $ — > —n and 
the excited boundary state to $ — ► n as x — > — 00. If however 7r/2 < <fo < ^ there is no 
boundary excited state: $ is too near either of the bulk vacuua values ±n at the boundary 
to drift over to the other as x — > — 00. 

In the fixed degenerate case (i.e. <po = 0), the interpretation is much the same. <po = 
is exactly between the bulk ground state field configurations $ = ±7r. Thus $ — > ±7r as 
x — > —00 are energetically equivalent. Because of this degeneracy there is the possibility 
that R^a does not vanish if a 7^ (3. To derive these cross-scattering matrix elements, one 
would presumably solve 

f(M ± )B a = ^R b a i(9)Al(-9)B p , (4.5) 

where / is a combination of fields appearing in ( p.!2|) . Thus the matrix elements diagonal 
in a and j3 remain unchanged. With the presence of scatttering between boundary states, 
the unitarity condition that must be satisfied is 

R b a i(e)RtJ(-6) = 8 ac 8 aj . (4.6) 

Given that we know in our case (as is easily checked) the diagonal matrix elements satisfy 
unitarity by themselves, it is straightforward to check that unitarity forces R%% oc 8 a p. 
Thus we still have no scattering between boundary states. The construction of the two 
boundary states B a in this case will include a zero momentum soliton (if the pole appears 
in P_|_) or anti-soliton (if the pole appears in P-). 

In the free (degenerate) case the interpretation of the boundary state differs. Con- 
sidering the massless limit, the twofold degenerate ground states are characterized by 
asymptotic behaviour $ — > ±7r as t —>■ 00 and $ =F7r as t — > —00. These two regions of 
different regions of $ are separated by a domain wall which can be thought of as a zero 
momentum particle present in the boundary state. This description is analogous to the 
description of semi-infinite Ising model at criticality in [Q. 

Our discussion has not yet exhausted all the poles in the physical strip < 9 < in. 
For every pole in the range < 6 < in/2, there is a corresponding pole in the range 
in/2 < 9 < in. These additional poles correspond to scattering in the cross channel. 
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5. Matching Ghoshal and Zamolodchikov's Solution 

We can now match our boundary S-matrix to that of Ghoshal and Zamolodchikov 
[IJ, thereby relating the formal parameters in the latter to the physical parameters in the 
action ( 2.1CQ . At the free-fermion point (corresponding to the limit A — > 1 in the notation 
of 0), the solution of the boundary Yang-Baxter equation yields 

p+(0) = cos(£ - id)R(9); 

p-(0) = cos(£ + id)R(d); 

Q+(9) = -i*± S mh(29)R(9); ^ 

Q-(e) = -z— sinh(26>)i?(0); 

where k± and £ are free parameters. Using the aforementioned gauge transformation, we 
set k + = k- = k. From the results of the last section we thus see that the phase a removed 
from k± in 0] is nothing other than the phase which specifies the conformal free boundary 
condition. 

The function R{6) is determined through the boundary unitarity and boundary cross- 
unitarity equations. In the free-fermion limit these constraints yield 

R{9) = ^7<t{ti, -i0)<r(W, ~ie), (5.2) 
where the parameters r\ and •& are related to k and £ through the equations 



cos(?7) cosh(#) = — — cos(£); 

k 

cos 2 (?7) + cosh 2 (i?) = 1 + — . 



(5.3) 



The function a(x,u) is found to be II 

a(x, u) = cosx (2 cos(tt/4 + x/2 - u/2) cos(tt/4 - x/2 - u/2))~ 1 . (5.4) 

With some algebra one can show that this solution of Ghoshal and Zamolodchikov 
maps onto our derivation of the boundary scattering matrix with the following identifica- 
tion of the free parameters: 

k= (l- 7 cos(^)+7 2 /4)" 1/2 ; 

cos£ = k (l - |cos</> ) ; (5.5) 

7 

sin£ = — k— sin</> . 
2 



5 The denominator of (5.23) in [0 should be Il 2 (a;, 7r/2)Il 2 (-x, -vr/2), and (5.24) should read 
a(x, u)a{x, —u) = cos 2 x[cos{x + \u) cos(x — Aw)] -1 . 
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In the limit 7 — > 0, we see that k = 1 and £ = yielding in Ghoshal and Zamolodchikov's 
solution, -P + (6>) = P~(9) and <3 + (6>) = Q~(9), in accordance with charge conjugation 
symmetry. In the limit of fixed boundary condition, 7 — > 00 we see that k = and 



forcing their amplitudes Q-t to vanish in accordance with U(l) charge conservation. We 
point out that the relationship between £ and <po is not in exact accordance with that 
conjectured in (see 5.27). There, it was taken that when </>o vanished, £ was forced to 
vanish in order to ensure P + (9) = P~{9). But £ = 7r also sets -P + (#) = P~(9). 

It should be emphasized that ( |5.6| ) depends strongly on the conventions used in writing 
the bulk Dirac action. Taking 7^ — > —7^ (where 7^ are the gamma matrices used in the 
action) leads to the relationship as conjectured in [0, namely 



In this case A = m and the poles in R(9) shift to reflect the ground states in the bulk are 
now $ = 2mr. 

For both ( |5.6| ) and ( |5.7| ) , taking $ — > — $ implements charge conjugation. However 
it is easy enough to construct Dirac actions where this action does not implement charge 
conjugation exactly, for example where the bosonisation relations are ip± = e ±lc?!> , ^ ± = 
_l_^ e T*0_ j n thjg case relation between </>o and £ takes the form <po = £±7r/2. Conventions 
of this sort are used in 0] . However these authors have taken the relation between £ and <po 
as given in ( |5.7| ) . This has led them to find a 'boundary anomaly'. This anomaly should 
vanish when the correct relation corresponding to their conventions between cpo and £ is 
used. 
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f = 4> + (2n + 1)tt, 



(5.6) 



£ = 4> + 2nir. 



(5.7) 
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